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Me-bhod for Perf onoing a Domain Transf orma-blon of a Dlgx-bal 
Signal from ^he Time Domain in'bo the Frequency Domain and 
Vice Versa. 

Cross Reference to Related Applications 

This application claims the benefit of priority of U.S. 
Provisional Application no 60/507,210, filed 29 September 
2003, and U.S. Provisional Application no 60/507,440, filed 
29 September 2003, the contents of each being hereby 
incorporated by reference in its entirety for all purposes. 

Further, the following commonly-owned applications are 
concurrently-filed herewith, and herein incorporated in its 
entirety: 

^^Method for Transforming a Digital Signal form the Time 
Domain into the Frequency Domain and Vice Versa, " Atty. 
Docket no. P100444, and 

^'Process and Device for Determining a Transforming 
Element for a Given Transformation Function, Method and 
Device for Transforming a Digital Signal from the Time Domain 
into the Frequency Domain and vice versa and Computer 
Readable Medium," Atty. Docket No. P100452. 

Background 

25 This invention relates generally to digital signal 

processing, and in particular to systems and methods for 
performing domain transformations of digitally-formatted 
signals . 

30 Domain transformations, for example the discrete cosine 

transform (DCT) , are widely used in modern signal processing 
industry. Recently, a variant of the DCT, called integer DCT, 
has attracted a lot of research interests because of its 
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important role in lossless coding applications. The term 
"lossless" means that the decoder can generate an exact copy 
of the source signal from the encoded bit-stream. 

The DCT is a real-valued block transform. Even if the input 
block consists only of integers,, the output block of the DCT 
can comprise non-integer components. For convenience, the 
input block is referred to as input vector and the output 
block as output vector. If a vector comprises only integer 
components, it is called an integer vector. In contrast to 
the DCT, the integer DCT generates an integer output vector 
from an integer input vector. For the same integer input 
vector, the integer output vector of integer DCT closely 
approximates the real output vector of DCT. Thus the integer 
DCT keeps all the good properties of the DCT in spectrum 
analysis . 

An important property of the integer DCT is reversibility. 
Reversibility means that there exists an integer inverse DCT 
20 (I DCT) so that if the integer DCT generates an output vector 
^ from an input vector x, the integer IDCT can recover the 
vector X from the vector Sometimes the integer DCT is also 
referred to as the forward transform, and the integer IDCT as 
the backward or inverse transform. 

25 

A transform called integer modified discrete cosine transform 
(IntMDCT) is recently proposed and used in the ISO/IEC MPEG-4 
audio con^pression. The IntMDCT can be derived from its 
prototype - the modified discrete cosine transform (MDCT) . 
30 The disclosure by H. S. Malvar in ^^Signal Processing with 

Lapped Transforms'' Artech House, 1992 provides an efficient 
realization of MDCT by cascading a bank of Givens rotations 
with a DCT-IV block. It is well known that Givens rotation 
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can be factorised into three lifting steps for mapping 
integers to integers. See e.g., R. Geiger, T. Sporer, J. 
Keller, K. Brandenburg, ^^Audio Coding based on Integer 
Transforms" AES 111^^ Convention, New York, USA, Sept. 2001. 
5 Integer transforms can be directly converted from their 
prototypes by replacing each Givens rotation with three 
lifting steps. Because in each lifting step there is one 
rounding operation, the total rounding number of an integer 
transfom is three times the Givens rotation number of the 

10 prototype transform. For discrete trigonometric transforms 
(for example the Discrete Fourier Transform (DFT) or the 
Discrete Cosine Transform (DCT) ) , the number of Givens 
rotations involved is usually at Nlog2N level, where N is 
the size of the blocks, i.e. the amount of data symbols 

15 included in each block, the digital signal is divided into. 
Accordingly, the total rounding number is also at JVlogjJV 
level for the family of directly converted integer 
transforms. Because of the roundings, an integer transform 
only approximates its floating-point prototype. The 

20 approximation error increases with the number of roundings. 



What is therefore needed is a method and system 
implementation for performing domain transformations using 
fewer rounding operations, thereby resulting in more accurate 
25 and less computationally intensive transformations. 

Summary of the Inven-bion 

The invention solves the problem of providing a method for 
performing a domain transformation of a digital signal from 
30 the time domain into the frequency domain and vice versa 

which comprises a significantly reduced number of rounding 
operations. 



wo 2005/031595 PCT/SG2004/000120 

4 

In one embodiment of the invention, a method for performing a 
domain transformation of a digital signal from the time 
domain into the frequency domain and vice versa includes 
performing the transformation by a transforming element which 
5 includes a plurality of lifting stages, wherein the 

transformation corresponds to a transformation matrix. At 
least one lifting stage includes at least one auxiliary 
transformation matrix, which is the transformation matrix 
itself or the respective transformation matrix of lower 
10 dimension. Further, each lifting stage includes a rounding 
unit. The method further includes processing the signal by 
the respective rounding unit after the transformation by the 
respective auxiliary transformation matrix. 

15 Brief Descrip-bion of -the Drawings 

Figure 1 shows the architecture of an audio encoder 
according to an embodiment of the invention; 

Figure 2 shows the architecture of an audio decoder 
according to an embodiment of the invention, which 
20 corresponds to the audio coder shown in figure 1; 

Figure 3 shows a flow chart of an embodiment of the 
method according to the invention; 

Figure 4 illustrates an embodiment of the method 
according to the invention using DCT-IV as the transformation 
25 function; 

Figure 5 illustrates the algorithm for the reverse 
transformation according to the embodiment of the method 
according to the invention illustrated in figure 4; 

Figure 6 shows a flow chart of an embodiment of the 
30 method according to the invention; 

Figure 7 illustrates an embodiment of the method 
according to the invention using DCT-IV as the transformation 
function; 
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Figure 8 illustrates the algorithm for the reverse 
transformation according to the embodiment of the method 
according to the invention illustrated in figure 7; 

Figure 9 shows the architecture for an image archiving 
5 system according to an embodiment of the invention; 

Figure 10 illustrates an embodiment of the method 
according to the invention using DWT-IV as the transformation 
function; 

Figure 11 illustrates the algorithm for the reverse 
10 transformation according to the embodiment of the method 
according to the invention illustrated in figure 10; 

Figure 12 illustrates a method for performing a domain 
transformation of a digital signal in accordance with one 
embodiment of the invention; 
•15 Figure 13 shows forward and reverse transform coders 

used to assess the transform accuracy of presented DCT and 
FFT transformation methods in accordance with the present 
invention . 

20 Detailed Description of Specific Embodiments 

Figure 12 illustrates a method for performing a domain 
transformation of a digital signal in accordance with one 
embodiment of the present invention- Initially at 1210, a 
transformation of the digital signal is performed by a 

25 transforming element comprising a plurality of lifting 

stages. The transformation corresponds to a transformation 
matrix wherein at least one of the lifting stages comprises 
an auxiliary transformation matrix and a rounding unit, the 
auxiliary matrix comprising the transformation matrix itself 

30 of the corresponding transformation matrix of lower 

dimension. Exemplary embodiments of the process 1210 and 
domain transformation which can be used in the domain 
transformation are further described below. 
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Next at 1220, a rounding operation of the signal is performed 
by the rounding unit after the transformation by the 
auxiliary matrix. 

According to preferred embodiments of the method according to 
the invention, the data symbols of the digital signal are 
provided to the transforming element as a data vector. In 
each lifting stage, the data vector or a part of the data 
vector is transformed by a domain transformer, and the result 
of the transformation is rounded subsequently to an integer 
vector. This is in contrast to any method according to the 
state of the art, wherein the rounding process is performed 
within the transformation process for the individual element 
or data symbol in each block of the digital signal. Thus, the 
number of rounding operations in the method according to the 
invention is greatly reduced. Due to the reduced number of 
rounding operations, the method according to the invention 
does not require large computation time and computer 
resource- Further, the approximation error of integer domain 
transformations can be significantly reduced. 

In a preferred embodiment, the invention provides a method 
for realizing the integer type-IV DCT transformation. The 
method according to the invention requires a significantly 
reduced number of rounding operations compared to prior art 
methods. As a result, the approximation error is greatly 
reduced, in the case of DCT- IV it can be reduced from the 
usual JVlogoiV level to be as low as 1.5N in one embodiment 
and as low as 2.5N in another embodiment, where N denotes the 
block size of the digital signal. The method according to the 
invention is low in computatdcnrval complexity and modular in 
structure . 
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The method and the device according to the invention can be 
used for any types of digital signals, such as audio, image 
or video signals. The digital signal, which is a digitised 
5 signal, corresponds to a physical measured signal, may be 

generated by scanning at least one characteristic feature of 
a corresponding analog signal (for example, the luminance and 
chrominance values of a video signal, the amplitude of an 
analog sound signal, or the analog sensing signal from a 
10 sensor) . The digital signal comprises a plurality of data 

symbols. The data symbols of the digital signal are grouped 
into blocks, with each block having the same predefined 
number of data symbols based on the sampling rate of the 
corresponding analog signal. 

15 

The method according to the invention can be used for 
transforming an input digital signal which represents integer 
values to an output signal which also represents integer 
values. The transformation method according to the invention 

20 is reversible. The output signal may be transformed back to 
the original input signal by performing the transformation 
method according to the invention. Such a reversibility 
property of the transformation according to the method of the 
invention can be used in lossless coding in which the output 

25 signal should be identical to the original input signal. 

Such an integer transformation of signals according to the 
invention can be used in many applications and systems such 
as MPEG audio, image and video compression, JPEG2000 or 
30 spectral analyzers (for analyzing Infrared, Ultra-violet or 
Nuclear Magnetic Radiation signals) . It can also be easily 
implemented in hardware systems such as in a fixed-point 
Digital Signal Processox (DSP) , without having to consider 
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factors such as overflow in the case of a real-value signal 
transformation . 

A preferred embodiment of the invention described below is 
5 suitable for both mono and stereo applications. In mono 

applications, two consecutive blocks of samples are groupeci 
and processed together. This introduces signal delay of one 
block length when compared to single-block processing. 
However, in stereo applications, this extra block delay can 
10 be prevented, if simultaneous sample blocks from the left and 
the right channel are grouped and processed together. 

According to the method according to the invention, the 
digital signal is transformed to the frequency domain by a 
15 transforming element. 

Preferably, the transforming element comprises a plurality of 
lifting stages. 

20 The transforming element can be illustrated based on the 

model of a lifting ladder. The lifting ladder model has two 
side pieces, each for receiving one of two groups of data 
symbols. Two or more cascading lifting stages are provided 
between the two side pieces. Each lifting stage receives a 

25 signal at one end (input end), and outputs a signal at the 

other end (output end) via a summation unit. A rounding unit 
is arranged at the output end. The lifting stages are 
arranged between the side pieces in an alternating manner, 
such that the output (or input) ends of adjacent lifting 

30 stages are connected to the different side pieces. 



It should be noted that although the transforming element is 
described in the form of the lifting ladder model, ii: is only 



wo 2005/031595 



PCT/SG2004/000120 



9 



to illustrate the transformation paths of the transforming 
element. However/ the invention shall not be limit.ed to said 
ladder model. 

5 Discrete cosine transforms, discrete sine transfoorms, 

discrete Fourier transforms or discrete W transforms are 
examples of transformations that may be performed by the 
method according to the invention. 

10 In one aspect of the invention, each lifting stage 

corresponds to a lifting matrix which is a block-triangular 
matrix comprising four sub-matrices, with two invertible 
integer matrices in one diagonal; An example of a lifting 
matrix is : 



wherein is the transformation function matrix^ in this 
case the DCT-IV transformation matrix, and Pl^ and P2j^ are 

20 permutation matrices. A peinmutation matrix is a matrix which 
changes the position of the elements in another matrix. The 
permutation matrices may also be identical. The invertible 
integer matrices may also be arranged in the lifting matrix 
in another way, for example as 



15 




(1) 



25 



• 0 
P2 




(2) 



30 



Preferably, the invertible integer matrices are diagonal 
matrices which comprise only components which are one or 
minus one. 
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Figure 1 shows the architecture of an audio encoder 100 
according to an embodiment of the invention. Thie audio 
encoder 100 comprises a conventional perceptual base layer 
5 coder based on the modified discrete cosine transform (MDCT) 
and a lossless enhancement coder based on the integer 
modified discrete cosine transform (IntMDCT) • 

An audio signal 109 which^ for instance, is provided by a 
10 microphone 110 and which is digitalized by a anaiog-to- 

digital converter 111 is provided to the audio encoder 100. 
The audio signal 109 comprises a plurality of da-ta symbols. 
The audio signal 109 is divided into a plurality of blocks, 
wherein each block comprises a plurality of data symbols of 
15 the digital signal and each block is transformed by a 

modified discrete cosine transform (MDCT) device 101. The 
MDCT coefficients are quantized by a quantizer 103 with the 
help of a perceptual model 102. The perceptual model controls 
the quantizer 103 in such a way that the audible distortions 
20 resulting from the quantization error are low. Thie quantized 
MDCT coefficients are • subsequently encoded by a bitstream 
encoder 104 which produces the lossy perceptually^ coded 
output bit stream 112. 

25 The bitstream encoder 104 losslessly compresses its input to 
produce an output which has a lower average bit-rate than its 
input by standard methods such as Huffman-Coding or Run- 
Length-Coding. The input audio signal 109 is also fed into an 
IntMDCT device 105 which produces IntMDCT coefficients. The 

30 quantized MDCT coefficients ^ which are the output of the 

quantizer 103, are used to predict the IntMDCT coefficients. 
-The quantized MDCT coefficients are fed into an inverse 
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quantizer 106 and the output (restored or non-quantized MDCT 

coefficients) is fed into a rounding unit 107. 

The rounding unit rounds to an integer value MDCT 
coefficients, and the residual IntMDCT coefficients, which 
are the difference between the integer value MDCT and the 
IntMDCT coefficients, are entropy coded by an entropy coder 
108. The entropy encoder, analogous to the bitstream encoder 
104, losslessly reduces the average bit-rate of its input and 
produces a lossless enhancement bitstream 113. The lossless 
enhancement bit stream 113 together with the perceptually 
coded bitstream 112, carries the necessary information to 
reconstruct the input audio signal 109 with minimal error. 

Figure 2 shows the architecture of an audio decoder 2 00 
comprising an embodiment of the invention, which corresponds 
to the audio coder 100 shown in figure 1. The perceptually 
coded bitstream 207 is supplied to a bitstream decoder 201, 
which performs the inverse operations to the operations of 
the bitstream encoder 104 of Figure 1, producing a decoded 
bitstream- The decoded bitstream is supplied to an inverse 
quantizer 202, the output of which (restored MDCT 
coefficients) is supplied to the inverse MDCT device 203. 
Thus, the reconstructed perceptually coded audio signal 209 
is obtained. 

The lossless enhancement bitstream 208 is supplied to an 
entropy decoder 204, which performs the inverse operations to 
the operations of the entropy encoder 108 of figure 1 and 
30 which produces the corresponding residual IntMDCT 

coefficients. The output of the inverse quantizer 202 is 
rounded by a rounding device 205 to produce integer value 
MDCT coefficients. The integer value MDCT coefficients are 
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added to the residual IntMDCT coefficients, thus producing 
the IntMDCT coefficients. Finally, the inverse IntMDCT is 
applied to the IntMDCT coefficients by an invexse IntMDCT 
device 206 to produce the reconstructed losslessly coded 
5 audio signal 210. 

Figure 3 shows a flow chart 300 of an embodiment of the 
method according to the invention using the DCT-IV as a 
transformation and using three lifting stages, a first 

10 lifting stage 301, a second lifting stage 302 and a third 
lifting stage 303. This method is preferably used in the 
IntMDCT device 105 of figure 1 and the inverse IntMDCT device 
206 of figure 2 to implement IntMDCT and inverse IntMDCT, 
respectively. In Fig. 3, Xi and X2 are first and second blocks 

15 of the digital signal, respectively, is an intermediate 

signal, and j^i and ^2 output signals corresponding to the 

first and second block of the digital signal, orespectively . 



As explained above, the DCT-IV-algorithm plays an important 
20 role in lossless audio coding. 

The transformation function of the DCT-IV comprises the 
transformation matrix Cj^ . According to this embodiment of 

the invention, the transforming element corresponds to a 
25 block— diagonal matrix comprising two blocks, wtierein each 

bLo.ck— comprises the transformation matrix Cjj^ . 



So, in this embodiment, the matrix corresponding to the 
transforming element is: 



30 



(3) 
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Cj^ shall in the context of this embodiment be referred to 
as the transformation matrix henceforth. 

5 The nvmiber of lifting matrices, and hence the number of 
lifting stages in the transformation element, in this 
embodiment of the invention, wherein DCT-IV is the 
transformation function, is three, 

10 The DCT-IV of an N-point real input sequence x(n) is defined 
as follows: 



y(m) 



[2^ / 



(m+i/2)(n+l/2)n 



N 



(4) 



15 Let Cf, be the transformation matrix of DCT-IV, that is 



(5) 



The following relation holds for the inverse DCT-IV matrix: 



20 



(6) 



With x^[x(n)l 

equation (4) can be 

25 expressed as 



y=cj:^x (7) 
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Now, let Xi/ ^2 be two integer iV'xl column vectors-. The 
column vectors Xir ^2 correspond to two blocks of the digital 
signal which, according to the invention, are transformed by 
one transforming element. The DCT-IV transforms of xj, x^ are 
5 jfi/ Y2f respectively. 



(8) 
(9) 



10 Combining (8) and (9) : 









A. 







(10) 



The above diagonal matrix is the block-diagonal matrix that 
15 the transforming element corresponds to- 



It is within the scope of the invention if the above equation 
is changed by simple algebraic modifications like the one 
leading to 




(11) 



Let be the counter diagonal matrix in (11), that is 



25 



(12) 



The matrix Tljy can be factorised as follows: 
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(13) 



where is the JSTxiV^ identity matrix- 

Equation (13) can be easily verified using the DCT— IV 
property in (7). Using (13), Equation (11) can be expressed 
as 







[ In 1 









































(14) 



The three lifting matrices in equation (14) correspond to the 
three lifting stages shown in Fig. 3. Each lifting matrix in 
(14) comprises the auxiliary matrix , which is the 
transformation matrix itself. 

From (14), the following integer DCT-IV algorithm tliat 
computes two integer DCT-IVs with one transforming element is 
derived . 

Figure 4 illustrates the embodiment of the method according 
to the invention using DCT-IV as the transf oinnation function. 
This embodiment is used in the audio coder 100 shown in Fig. 
1 for implementing IntMDCT. Like in Figure 3, Xi and are 
two blocks of the input digital signal, z is an intermediate 
signal, and and are corresponding blocks of th.e output 
signal. 

The three lifting stages illustrated in figure 4 correspond 
to the three lifting matrices in equation (14) . 
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As illustrated by figure A, the time to frequency domain 
integer transform is determined by the following: 

In the first stage 401, X2 is transformed by a DCT-IV 
transformation 402 and the DCT-IV coefficients are rounded 

403. The rounded DCT-IV coefficients are then added to xj 

404. Thus, the intermediate signal z is generated. So, the 
intermediate signal z fulfils the equation: 



In the second stage 405, z is transformed by a DCT-IV 
transformation 406 and the DCT-IV coefficients are rounded 
407. From the rounded DCT-IV coefficients xx is then 
subtracted. Thus, the output signal j^j is generated. So, the 
output signal jfi fulfils the equation: 



In the third stage 409, y.i is transformed by a DCT-IV 
transformation 410 and the DCT-IV coefficients are rounded 
411. The rounded DCT-IV coefficients are then subtracted from 
z. Thus, the output signal ^2 is generated. So, the output 
signal ^2 fulfils the equation: 



(15a) 



(15b) 



y2-Ac^yr\ 



(15c) 



where [*J denotes rounding operation. 



Figure 5 illustrates the algorithm for the reverse 
transformation according to an embodiment of the method 
according to the invention using DCT-IV as the transformation 
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function. This embodiment is used in the audio decoder 200 
shown in Fig. 2 for implementing inverse IntMDCT. The 
algorithm illustrated in figure 5 is the inverse of the 
algorithm illustrated in figure 4. The denotations for the 
different signals ^2, Xi, X2 and z are chosen 
corresponding to the denotations of figure 4. 

As illustrated by Figure 5, the frequency to time domain 
integer transform is determined by the following: 

In the first stage 501^ is transformed by a DCT— IV 
transformation 502 and the DCT-IV coefficients are rounded 

503. The rounded DCT-IV coefficients are then added to ^2 

504. Thus, the intermediate signal z is generated. So, the 
intermediate signal z fulfils the equation: 

^^feZlJ+Zi (16a) 

In the second stage 505, z is transformed by a DCT— IV 
transformation 506 and the DCT-IV coefficients are rounded . 
507. From' the rounded DCT-IV coefficients Zi is then 
subtracted. Thus, the signal X2 is generated. So, the signal 
X2 fulfils the equation: 




(16b) 



In the third stage 509, is transformed by a DCT-IV 
transformation 510 and the DCT-IV coefficients are rounded 
511. The rounded DCT-IV coefficients are then subtr-acted from 
z. Thus, the signal Xi is generated. So, the signal Xj 
fulfils the equation: 
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It can be seen that the algorithm according to the equations 
5 (16a) to (16c) is inverse to the algorithm according to the 
equations (15a) to (15c) . Thus, if used in the encoder and 
decoder illustrated in figures 1 and 2, the algorithms 
provide a method and an apparatus for lossless audio coding. 

10 The equations (15a) to {15c) and (16a) to (16c) further show 
that to compute two NxiN integer DCT-IVs, three NxN DCT- 
IVs, three Nxl roundings, and three Nxl additions are 
needed. Therefore, for one NxN integer DCT-IV, the average 
is: 

15 

i?C(JS0=l-5JNr (17) 
AC(N) - 1.5AC(C^ ) + 1.5JV^ (18) 

20 where RC(.) is the total rounding number, and AC(,) is the 
total number of arithmetic operations. Compared to the 
directly converted integer DCT-IV algorithms, the proposed 
integer DCT-IV algorithm reduces RC from level JVlogjiNT to N, 

25 As indicated by (18), the arithmetic complexity of the 

proposed integer DCT-IV algorithm is about 50 percent more 
than that of a DCT-IV algorithm. However, if RC is also 
considered, the combined complexity (AC+J^C) of the proposed 
■^ai-gorithmi doBS' not much exceed that of the directly converted 

30 integer algorithms. Exact analysis of the algorithm 
complexity .depends on the DCT-IV algorithm used. 
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As shown in figures 4 and 5, the proposed integer DCT-IV 
algorithm is simple and modular in structure. It can use any 
existing DCT-IV algorithms in its DCT-IV computation block. 
The proposed algorithm is suitable for applications that 
5 require IntMDCT, e.g. in the MPEG-4 audio extension 3 
reference model 0. 

Figure 6 shows a flow chart 600 of an embodiment of the 
method according to the invention using five lifting stages^ 
a first lifting stage 601, a second lifting stage 602, a 
third lifting stage 603, a fourth lifting stage 604 and a 
fifth lifting stage 605. This method can be used in the 
IntMDCT device 105 of figure 1 and the inverse IntMDCT device 
206 of figure 2 to implement IntMDCT and inverse IntMDCT, 
respectively. In Fig. 6, Xi X2 are first and second blocks 

of the digital signal, respectively, zi, z^, z^ are 
intermediate signals, and y^x and Y2 are output signals 
corresponding to the first and second block of the digital 
signal , respectively . 

Figure 7 shows a flow chart of an embodiment of the method 
according to the invention, wherein the transformation 
function is a DCT-IV transformation function. 

25 The DCT-IV matrix (see equation (5)) can be factorised 

as 

C^^R^P^ (20) 

30 where is an even-odd matrix, i.e., a permutation inatrix 
reordering the components of the vector 
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X ^ 



x{0) 
x(N-l) 



by separating the components which correspond to even indices 
from the components which correspond to odd indices, i.e., 
5 such that 



^4-1) 

^2(0) 



(21) 



10 



where further 



1 = 2^2^ 



^1 In/2 



^2 



'N/2 



^3 ^N/2 



(22) 



with 



15 = "" (^Ar/2-Pjy/2 '^' ^^n/2]Pn/2 



^NI2 



(23) 
(24) 
(25) 



20 



and 



^N/2 




^ N/2 ^2 




^N/2 


^1 ^Nl.2. 




^N/2 




^2 ^ N/2 



(26) 
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-tan 



- tan- 



- tan- 



L 



3n 



SN 



(27) 



and 



sin 

4N 



7t 



4N 



3n 
4N 



(28) 



Equation (20) can be further written as 



10 CjJ ^R^R2ST^T^P^ 



(29) 



where 



15 



20 



(30) 



The transforming element used in the embodiment illustrated 
in figure 7 corresponds to the equation (29) . Each lifting 
matrix Sr T2 and T3 in (29) comprises the auxiliary matrix 



which is the transformation matrix itself. 



The transforming element compri'siers" five lifting stages which 
correspond to the five lifting matrices of equation (29) . 
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Further, the transforming element comprises a data shuffling 
stage corresponding to the permutation matrix . 

In figure 7, the input of the first lifting stage are the two 
5 blocks of the digital signal xi and X2, zx^ Z2 and £3 are 
intermediate signals and j/i and jf2 are output signals 
corresponding to the first and second block of the digital 
signal , respectively . 

10 The input of the transforming element x and the two input 

blocks of the first lifting stage of the transforming elennent 
Xi and K2 fulfil the equation 




(31) 



15 

In the following, the first lifting stage 701 is explainecl, 
which is the lifting stage corresponding to the lifting 
matrix . 

be the output vector of the first lifting stage, for 
the time being without rounding to integer values, i.e. 




(32) 



25 Using the definition of provided by equation (22), 
equation (29) is re-written as 



20 Let 
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2 






h.' 




















2 







(33) 



where denotes the identity matrix of dimension N/2. 

2 

Since in this embodiment, a reversible algorithm for integer 
DCT-IV in this provided, there is included a rounding to 
integer values. So, according to equation (30), in the first: 
step 706 of the first lifting stage 701, xx is multiplied hy 
. The result of this multiplication is rounded to integer- 
values in step 707. The rounded values are then added to 
in step 708. Thus, the intermediate signal zi fulfils the 
equation 




(34) 



where [*J denotes rounding operation. 

Since the other four lifting stages 702, 703, 704, 705 of the 
transforming element illustrated in figure 7, which 
correspond to the matrices I^, S, and respectively, 

have identical structures as the first lifting stage 701, 
their explanation is omitted. It only has to be noted that in 
the adding step 709 of the second lifting step 702, Xi is, 
according to the definition of , multiplied with -D^,^ . 



In the following, the lifting stages of the transforming 
element of the reverse transformation are explained with 
reference to figure 8. 
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Figure 8 illustrates the lifting stages of the transforming 



element for the reverse transformation of the transformation 
illustrated in figure 1. 

5 In figure 8, the input of the first lifting stage are the two 
blocks of the digital signal and y2, ^i/ £2 and Z3 are 
intermediate signals and xj and are output signals 
corresponding to the first and second block of the digital 
signal, respectively . 



The last lifting stage 805 illustrated in figure 8 is inverse 
to the first lifting stage 701 illustrated in figure 7. So, 
in the first step 806 of the last lifting stage 805, ±s 
multiplied by . The result of this multiplication is 

15 rounded to integer values in step 807 • The rounded values are 
then subtracted from zx in step 808. Thus, the signal 
fulfils the equation 



Since the other four lifting stages 801, 802, 803, 804 of the 
transforming element illustrated in figure 8, which are 

25 inverse to the lifting stages 705, 704, 703 and 702 of the 
t-ranrsf orming element illustrated in figure 7 respectively, 
have identical structures as the last lifting stage 805, 
their explanation is omitted. It only has to be noted that 
after— the— adding step 809 of the fourth lifting step 804^ the 

30 result of the adding step 809 is multiplied with ~jD^/2 to 

produce Xi. 



10 




(35) 



20 



where [*J denotes rounding operation. 



wo 2005/031595 PCT/SG2004/000120 

25 

It can be seen that the lifting stages 805, 804, 803, 802 and 
801 of figure 8 are inverse to the lifting stages 701 to 705 
of figure 1 , respectively. Since also the permutation of the 
input signal corresponding to the matrix can be reversed 

5 and an according data shuffling stage is comprised by the 
inverse transforming element, the provided method is 
reversible, thus, if used in the audio encoder 100 and the 
audio decoder 200 illustrated in figures 1 and 2, provi^ding a 
method and an apparatus for lossless audio coding. 

10 

An analysis of the number of roundings used in this 
embodiment is given at the end of the description of the 
invention. 

15 Figure 9 shows the architecture for an image archiving system 
according to an emboddLment of the invention. 

In figure 9 an image source 901, for instance a camera, 
provides an analog image signal. The image signal is 

20 processed by a analog-to-digital converter 902 to provicie a 
corresponding digital image signal. The digital image signal 
is losslessy encoded by a lossless image encoder 903 whdLch 
includes a transformation from the time domain to the 
frequency domain. In this embodiment, the time domain 

25 corresponds to the coordinate space of the image. The 

lossless coded image signal is stored in a storage device 
904, for example a hard disk or a DVD. When the image is 
needed, the losslessly coded image signal is fetched from the 
storage device 904 and provided to a lossless image decoder 

30 905 which decodes the losslessly coded image signal and 

reconstructs the original image signal without any data loss. 
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Such lossless archiving of image signals is important, for 
example, in the case that the images are error maps of 
semiconductor wafers and have to be stored for later 
analysis. 

In the following, a further embodiment of the method for the 
transformation of a digital signal from the time domain to 
the frequency .domain and vice versa according to the 
invention is explained. This embodiment is preferably used in 
the lossless image encoder 903 and the lossless image decoder 
905 of the image archiving system illustrated in figure 9. 

Figure 10 illustrates an embodiment of the method according 
to the invention using DWT~IV as the transformation function. 

The DWT-IV of a W--point real input sequence x(n) is defined 
as follows: 



, . [2^ , ^ . (k (m + l/2)(n+l/2)27r\ ^ ^ 



(36) 



The transformation matrix Wj^ of DWT-IV is given by 



[m + l/2)(n + y2)2n ' 
N 



(37) 



The DWT-IV matrix is factorised into the following form: 



(38) 



R„ is a N-kN rotation matrix defined as 
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1 r / 



1^/2 ^ N/2 



N/2 ^N/2 



(39) 



^N/2 is the identity matrix of order N/2. J^^^ is the counter 
identity matrix of order N/2, i.e.. 



N^ 
2 



(40) 



Pjy is a NxN permutation matrix 



Pn 



^N/2 



N/2 



(41) 



T is a NxN matrix, given by 



1 

V2 



C" D 



D 



(42) 



where Cj^/a is the DCT-IV matrix of order N/2, i.e., a 
corresponding transformation matrix of lower dimension. 



/ (m + V2)(n + V2)7r y 



(43) 



is an order N/2 diagonal matrix given by 
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-1 



D 



N/2 



(44) 



-1 



Rj, and T can further be factorised into a product of lifting 
matrices : 



^Nt2 








^NI2 


^1 




^Ni2 




^3 ^N/2 



(45) 



where = (>/2/^-C^^^D^)c^ ^ = 



-C 



J2 



JV 

^3 ^-^^iV/2-^JV/2 ""'';sr/2 / 



10 



and 



^N/2 










^N/2 


El 


^N/2 

■J 




^N/2 




^3 ^N/2 



(46) 



where 



15 Hj^^H^'^ -tan(;r/8)- J 



jy/2 



-0.414 



-0.414 



-0.414 



(47) 



and 



^ = sm(7r/4)-7j^ = 



0.707 



0.707 



0.707 



20 



(48) 



Thus, equation (38) can be written in the forrm 
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(49) 



The lifting matrices ^ and i; can be aggregated to the 
5 lifting matrix S : 



(50) 



From (49) and (50), the following factorisation formula for 
10 the DWT-IV matrix can be obtained: 



C^=^ R1R2S. ^2^3-^ 



(51) 



Equation (51) indicates that a transforming element for the 
15 integer DWT-IV transform can be used which comprises five 
lifting stages. Each lifting matrix S, T2 and T3 in (51) 
comprises the auxiliary matrix Cj^^2 f which ±s the 
transformation matrix itself. 



20 Further, the transforming element comprises a data shuffling 
stage corresponding to the permutation matrix - The data 
shuffling stage rearranges the components order in each input 
data block. According to P^^r the input data vector is 
rearranged in the following way: the first half of vector 

25 remains unchanged; the second half of vector is flipped, 
i.e. : 
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1 
















^JV/2 



























(52) 



In figure 10, the input of the first lifting stage are the 
two blocks of the digital signal Xj and X2, £if and £3 are 
5 intermediate signals and and ^2 are output signals 

corresponding to the first and second block of the digital 
signal , respectively . 

The input of the transforming element x and the two input 
10 blocks of the first lifting stage of the transforming element 
Xi and X2 fulfil the equation 



(53) 



15 In the following, the first lifting stage 1001 is explained, 
which is the lifting stage corresponding to the lifting 
matrix . 



Let [— be the output vector of the first lifting stage, for 
L— J 

.20 the time being without rounding to integer values, i.e. 



X2 



(54) 
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Using the definition of provided by equation (45), 
equation <51) is re-written as 



^2 



^2 



(55) 



Since in this embodiment/ a reversible algorithm for integer 
DWT-IV in this provided^ there is included a rounding to 
integer values. So, according to equation (55), in the first 
step 1006 of the first lifting stage 1001, Xj is multiplied 
10 by . The result of this multiplication is rounded to 

integer values in step 1007, The rounded values are then 
added to X2 in step 1008- Thus, the intermediate signal _2j 
fulfils the equation 

15 £i="|^'5,J+^ ^^^^ 
where [*J denotes rounding operation. 

Since the other four lifting stages 1002, 1003, 1004, 1005 of 
20 the transforming element illustrated in figure 10, which 

correspond to the matrices 7^, 5, R2 and respectively, 

have identical structures as the first lifting stage 1001, 
their explanation is omitted. It only has to be noted that in 
the adding step 1009 of the second lifting step 1002, Xj is, 
25 according to the definition of 2^, multiplied with ^^72- 

In the following, the lifting stages of the transforming 
element of the reverse transformation are explained with 
reference to figure 11. 

30 



10 



20 
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Figure 11 illustrates the lifting stages of the tr-ansf orming 
element for the reverse transformation of the trarxs formation 
illustrated in figure 10. 

In figure 11, the input of the first lifting stage are the 
two blocks of the digital signal yi and ^ir ^2 and Z3 are 
intermediate signals and xi and K2 are output signals 
corresponding to the first and second block of the digital 
signal , respectively . 



The last lifting stage 1105 illustrated in figure 11 is 
inverse to the first lifting stage 1001 illustrated in figure 
10. So, in the first step 1106 of the last lifting stage 
1105, Kx is multiplied by . The result of this 

15 multiplication is rounded to integer values in step 1107. The 
rounded values are then subtracted from zi in step 1108. 
Thus, the signal X2 fulfils the equation 



where [*J denotes rounding operation. 



Since the other four lifting stages 1101, 1102, 1103, 1104 of 
the transforming element illustrated in figure 11, which are 

25 inverse to the lifting stages 1005, 1004, 1003 and 1002 

respectively, have identical structures as the last lifting 
stage 1105, their explanation is omitted. It only Inas to be 
noted that after the adding step 1109 of the fourth lifting 
step 1104, the result of the adding step 1109 is mioltiplied 

30 with produce Xi- 
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It can be seen that the lifting stages 1105, 1104, 1103, 1102 
and 1101 of figure 11 are inverse to the lifting stages 1001 
to 1005 of figure 10, respectively. Since also the 
permutation of the input signal corresponding to the matrix 
csri be reversed and an according data shuffling stage is 

comprised by the inverse transforming element, the whole 
provided method is reversible, thus, if used in the lossless 
image encoder 903 and the lossless image decoder 90 5 
illustrated in figure 9, providing a method and an apparatus 
for lossless image coding. 

Although in the explained embodiments, the embodiments of the 
method according to the invention for DCT-IV were used for 
audio coding and the embodiment of the method according to 
the invention for DWT-IV was used for image coding, the 
embodiments of the method according to the invention for DCT- 
IV can as well be used for image coding, the embodiment of 
the method according to the invention for DWT-IV can as well 
be used audio coding and all can be used as well for the 
coding of other digital signals, such as video signals. 

Considering the equations (34) and (35) one sees that there 
are N/2 roundings in each lifting stage. Therefore, 
considering equation (26), one sees that the total rounding 
number for the transforming element for the embodiment of the 
DCT-IV algorithm illustrated in figures 7 and 8 according to 
the invention is five times N/2, that is 2.5Nf which is 
significantly lower than Nlog2N according to the state of the 
art. 

Again considering equation (29), one sees that the majority 
of computation power is used in the four N/2 point DCT-IV 
subroutines corresponding to a multiplication with Cj^,^ / when 



wo 2005/031595 PCT/SG2004/000120 

34 

N ±s a large value, e.g. N ^ 1024. Because the floating-point 
DCT-IV can be calculated using two half-length DCT-IV plus 
pre- and post-rotations, the arithmetic complexity of the 
proposed integer DCT-IV is roughly estimated to be twice that 
5 of the floating-point DCT-IV. 

A similar conclusion can be drawn for the integer DWr-IV 
transformation function. 



10 In the following, a further embodiment of the method for the 
transformation of a digital signal from the time domain to 
the frequency domain and vice versa according to the 
invention is explained. 

15 In this embodiment of the present invention, the domain 

transformation is a DCT transform, whereby the block size N 
is some integer and where the input vector comprises two sub- 
vectors . 



20 Let be the N kN transform matrix of DCT (also called 
Type-II DCT) : 

=V2;]v[^„cosMn+V2W^)] 



25 



where: 



VV2 if m-0 
1 if m^O 



(59) 



30 



and N is the transform size which can be a powex erf two, 
N = 2' in one embodiment and i is an integer greater than 
zero, m and n are matrix indices. 
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Let be the N^N transform matrix of type-IV DCT, as 
already defined above: 

Cir - VW[cos((/n + yiXn + V2>r/i\^)] ^ ^ 



10 



15 



As above, a plurality of lifting matrices will be used, ^hich 
lifting matrices are in this embodiment 2Nx2N matrices of 
the following form: 



^2N = 



(61) 



where is the xA^ identity matrix, is the NxN zexo 

matrix, Aj^ is an arbitrary NxN matrix. 



For each lifting matrix L,2n ^ ^ lifting stage reversible 
integer to integer mapping is realized in the same way as the 
2x2 lifting step described in the incorporated reference 
"Factoring Wavlet Transforms into Lifting Steps," Tech. 
20 Report, I- Daubechies and W. Sweldens, Bell Laboratories^ 
Lucent Technologies, 1996. The only difference is that 
rounding is applied to a vector instead of a single variable. 



In the above description of the other embodiments, it was 
25 already detailed how a lifting stage is realized for a 

lifting matrix, so the explanation of the lifting stages 
corresponding to the lifting matrices will be omitted in the 
following. 
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One sees that the transposition of Lj^^ , is also a 

lifting matrix. 



In this embodiment, the transforming element corresponds to a 
5 matrix, which is defined as a 2Nx2N matrix in the 

following way: 



o 



(62) 



The decomposition of the matrix T2jy into lifting matrices 
10 has the following form: 



-P3-L8-L7-L6-P2-L5-L4-L3-L2 L1P1 



(63) 



The matrices constituting the right hand side of the above 
15 equation will be explained in the following. 

PI is a first permutation matrix given by the equation 



20 



Pt: 



where J^^ is the NxN counter index matrix given by 



(64) 



0 


0 


0 


1 


0 


0 


1 


0 


0 




0 


0 


1 


0 


0 


0 



(65) 



Z5 and.-D^~is a NxN diagonal matrix v/ith diagonal element 
being 1 and -1 alternatively: 
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1 


0 


0 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 




0 


0 


0 


0 


0 


-1 



(66) 



P2 is a second peormutation matrix, an example of whd_ch is 
generated by the following MATLAB script : 

5 



Pd = eye(2*N) ; 
for i = 2:2:N, 

Pd(i,i) = 0; Pd(N+i,N+i) = 0; 
10 Pd(i,N+i) = 1; Pd{N+i,i) = 1; 

end 

Peo = zeros (2*N); 
for i = 1:N, 

Peo(i,2*i-l) = 1; 
15 Peo(i+N,2*i) = 1; 

end 

P2 = (Pd*Peo) ' ; 



20 As an example, when JV is 4, P2 is a 8x8 matrix given as 



1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


1 


0 


0 


0 


0 
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P3 is a third permutation matrix, an example of which is 
generated by the following MATLAB script: 



P3 = zeros (2*N) ; 
for i = 1:N, 

P3(i,2*i-1) = 1; 

P3 (N2-i+l,2*i) = 1; 

end 



10 



As an excimple, when N is 4, P3 is a 8x8matrix given as 



P3. 



1 


0 


0 


0 


0 


0 


0 


0- 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 



15 



LI is a first lifting matrix 



Ll = 



Zl 



■N 



for iSr-4 



(68) 



O 



AT 



(69) 



where 'Zlj^ is a Ni<N counter diagonal matrix given as: 



20 



Zl 



0 
0 
0 



0 
0 



— taii((2iv - i)n:/8N) 0 



0 

tan(3;r/8JV) 
0 
0 



- tan(;r/8JV)' 
0 
0 
0 



(70) 



L2 is a second lifting matrix: 
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L2 = 



(71) 



where Z2j, is a NxN counter diagonal matrix given as; 



Z2 



N 



0 0 0 sm((2Ar-l);r/4JSr)- 

0 0 ••■ 0 

0 sin{37c/4N) 0 0 

sin(;r/4iy^) 0 0 0 



(72) 



L3 is a third lifting matrix; 



10 



L3. 



I„ O 



■N 

Z3 



(73) 



where : 



Z3^=V2C;j' +Ijv+Zl^ 



15 



L4 is a fourth lifting matrix: 



(74) 



L4. 



o 



(75) 



20 where: 



Z4^=Cjr/>^ 



(76) 



L5 is a fifth lifting matrix: 



25 



L5 



(77) 
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where : 



Z5^--(>/2Cjr+I^) 



(78) 



L6 is a sixth lifting matrix: 



L6 = 



1.. 

Z6^ Ijf 



(79) 



10 where Z6f, is a NxN counter diagonal matrix given as; 



Z6 



0 0 0 tan(;r/8) 

0 0 tan(w/8) 0 

0 .* 0 0 

tan(;r/8) 0 0 0 



(80) 



L7 is a seventh lifting matrix: 



15 



L7 = 



Iat Z7^ 



(80.) 



where Z7„ is a NxN counter diagonal matrix given as: 



20 



Z7 



N 





0 


0 


0 


-sin(7r/4)' 




0 


0 




0 




0 


-sm(;r/4) 


0 


0 




-sm(;r/4) 


0 


0 


0 



(82) 



LB is an eighth lifting matrix: 



L8 = L6 



(83) 



5 



15 
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thus, resulting in the factorization as shown in (x) : 



«P3-L8 L7-L6-P2-L5-L4-L3 L2-L1-P1 (84) 



where Fl, V2 , and P3 are three permutation matrices. L; , / 
from 1 to 8, are eight lifting matrices. 

The lifting matrices 1-3, L4 and L5 comprise an auxiliary 
10 transformation matrix, which is, in this case, the 

transformation matrix C^^ itself. 

From Eq. (84), it is possible to compute the integer DCT for 
two input signals of dimension iVxl. 



As Eq. (84) provides a lifting matrix factorization which 
describes the DCT-IV transformation domain, its lifting 
matrices can be used in the manner shown herein to compute 
the domain transformation of an applied input signal. 

The equation (84) can be derived in the following way. 



The following decomposition can be derived using the 
disclosure from Wang, Zhongde, "On Computing the Discrete 
25 Fourier and Cosine Transforms", IEEE Transactions on 

Acoustics, Speech and Signal Processing, Vol. ASSP-33, No. 4 
October 1985: 
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c-=(B,r.(p,r 



(85) 



P -T 



is known^ wherein ^^^otes the transformation matrix of 

the discrete sine transform of type 2, 



D J 



is a NxN permutation matrix given by 



10 P^ = 



' Af/2 



H/2 



(86) 



n 



cos 



41) 



cos 



3jt 
4N 



- sin 



3jc 
4N 



— sxn 



4N 



(n - l) ir (n - l) at 

cos sin 



- sin 



4N 

(n - l) JC 
4N 



cos 



4N 

(n - l) 3t 

4N 



sm 



4N 



sm 



3jt 
4M 



cos 



3n 
4N 



cos 



4N 



15 and— 
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>/2 



1 -1 
1 1 



43 



1 -1 
1 1 



V2 



Equation (85) can be combined with the equation 



5 ""I^po 



(87) 



wherein P^o is an even-odd permutation matrix, 



'N/2 



'N/2 



N/2 



10 Rpo equals Tn, 




After transposition equation (87) converts to 

15 



^N/2 



^N/2 



^N/2 

^^72 



^N/2 



(88) 
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The combination of (85) and (88) yields 



*3 ^2 ^2 



1 

\/2 



C 



TV 

JV 
N/2 



^N/2 



Ri Pi 



W/2 



^N/2 
^N/2 



(89) 



where 



P2=(P^of •(Pz,)'=(Px>-P^of 
P3=P^ 

^2 



from (89), equation (84) can be easily derived. 

In this embodiment, the computation of the domain 
transformation requires only 4N rounding operations, as w±ll 
now be explained: 

Let be the number of real additions, ^(*) be the numbor 

of real multiplications, and y{^) be the number of real 
roundings, respectively. For the proposed IntDCT algoritlnm, 
one -getS4 

a(lntDCT)=lW+3a(DCT-Iv) 
/i(lntnCT) = 9iSr + 3fi(pCT - IV) 
y(EntDCT)-8J\r 

The above results are for two blocks of data samples, because 
the proposed IntDCT algorithm processes them tO-getharL. Thtas 
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for one block of data samples, the numbers of calculations 
are halved, which are 

ai(liitDCT) = 5.5JV^ +l^a(DCT - IV) 
fi, (intDCr) = 4SN + l,5/i(DCT - IV) 

yi(lntDCr) = 4JV^ 

where a^, /i^ , and are the number of real additions, number 
of real multiplications, and number of real roundings, for 
one block of samples, respectively. 

For DCT-IV calculation, the FFT-based algorithm described in 
the incorporated reference "Signal Processing with lapped 
Transforms," S. Malvar, Norwood, MA. Artech House, 1992, 
pp. 199-201 can be used, for which 

a (DCT - IV) = 1.5JSriog2 N 
Ai(DCT - IV) = 0.5J\riog2 N-hN 

Consequently : 

ai(lntDCT)» 2.25J\riog2 N^5.5N 
Ail (intDCr) = 0.75N logj iST + 6JSr 

In the following, a further embodiment of the method for the 
transformation of a digital signal from the time domain to 
the frequency domain and vice versa according to the 
invention is explained. 

In this embodiment a discrete fast fourier transform (FFT) is 
used as the domain transformation. 

Let F be the NxN transform m.atrix of the normalized FFT 
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(91) 



where ^ is the transform size which is some positive integer 
5 and, in one embodiment, can be a power of two, N ^2^ and * is 
an integer greater than zero. ^ and ^ are matrix indices. 

Under this embodiment, a permutation matrix P of dimension 
N is a matrix which includes indices 0 or 1. After 
10 multiplying it with a -^^1 vector (the matrix representation 
of the input signal) , the order of elements in the vector are 
changed . 

In this embodiment, lifting matrices are defined as 2N^x2N 
15 matrices of the following form: 



20 



Pi A 
O P, 



(92) 



where and ^2 are two permutation matrices, ^is the NxN 
zero matrix, A is an arbitrary matrix. For lifting 

matrix L, reversible integer to integer mapping is realized 
in the same way as the 2x2 lifting step in the 
aforementioned incorporated reference of I. Daubechies . As 
above, however, rounding is applied to a vector instead of a 
25 single variable. It is apparent that the transposition of 

is also a lifting matrix. 
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Further, let T be a ^NxlN transform matrix: 



O F 
F O 



(93) 



Accordingly, the modified transform matrix T (and 
accordingly the domain transformation itself) can be 
expressed as the lifting matrix factorization: 



T = 



10 



I 






F 




I O 


-QF 


°\ 


Co 


I 




F I 
. * J 



(94) 



where I is the identity matrix, and Q is a 

permutation matrix given as: 



15 



1 O 



O 



lxiV-1 
J 



( 95) 



and ^^-^ and ^*f-^ are row and column vectors of -^-1 zeros 
respectively . 

20 J is the counter index matrix given by 



(96) 
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In Eq, (96), blank space in the square bracket represents all 
zeros matrix elements. 

The lifting matrices on the right hand side of equation (94) 
comprise an auxiliary transformation matrix, which is, in 
this case, the transformation matrix F itself. 

A can be seen from Eq. (94), the lifting matrix factorization 
can be use to compute the integer FFT for two complex 
vectors using the methods as described herein. 

Under this embodiment, the computation of the domain 
transformation requires only 3N rounding operations, as will 
now be explained: 

Let:^(*) be the number of real additions, 

be the number of real multiplications, and 

be the number of real rounding operations, 
respectively. 

For the proposed IntFFT algorithm, we have 



a(liitFFT) = 6JN^ + 3a(FFT) 
^(lntFFT) = 3Ai(FFr) 



The~above results are for two 
the ..proposed IntFFT algorithm 
for one. block of data samples, 
are halved, which aire 

a, (IntFFT) « SiV^ + 1 .5a (FFT) 



blocks of data samples, because 
processes them together. Thus 
the nximbers of calculations 
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^(lntFFT)=1.5/i(FFr) 
yi(lntFFr)=3iV 

where , ^, and are the niimber of real additions^ nuunber 
of real multiplications, and niunber of real rounding 
5 opreations for one block of samples, respectively. 

For FFT calculation, the split-radix FFT (SRFFT) algoorithm 
can be used, for which: 

a(SRFFT) = SJVloga J\r -SJST + 4 

/x(SRFFT)=^^log2^^-3^^ + 4 



Consequently, we have: 

ai(liitFFT) = 4.5JV^log2 N -ISN + 6 
/^^(intFFT) -1.5JV^log2 N --A.SN + 6 

15 

Figure 13 shows forward and reverse transform coders used to 
assess the transform accuracy of the DCT transformation 
technique described above and the FFT domain transformation 
above. The test involved measuring the mean squared error 

20 (MSE) of the transform in accordance with the evaluation 

standards proposed by the MPEG-4 lossless audio coding group 
as described in "Coding of Moving Pictures and Audio: Work 
plan for Evaluation of Integer MDCT for FGS to Lossless 
Experimentation Framework," ISO/IEC JTC 1/SC 29/WG 11 N5578 

25 Pattaya, Thailand, Mar. 2003 incorporated herein. 

Specifically, the MSEs for IntDCT and integer inverse DCT 
(IntlDCT) are given as: 



30 



1 K-\ i N'\ 



(97) 
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where the error signal ^ is ^ for IntDCT, and ' for IntlDCT 

as in Fig 1. ^ is the total nimber of sample blocks used in 
the evaluation. 

5 The MSEs for IntFFT and integer inverse FFT (IntlFFT) are 
given as 



MSE = 




(98) 



10 where the error signal ^ is ^ for IntFFT, and ' for IntlFFT 

IHI 

as in Fig 1. " " represents norm of a complex value. ^ is the 
total number of sample blocks used in the evaluation. 

For both domain transformations, a total of 450 seconds with 
15 15 different types of music files are used in the 48 kHz/16- 
bit test set. Table I shows the test results. 

As can be seen from Table 1, the MSE generated using the 
systems and methods of the present invention is very minimal, 

20 and unlike conventional systems, is substantially independent 
of the processing block size. Referring to the DCT-IV domain 
transformation, the MSE only slightly increases with 
increasing block size N up to 4096 bits. The MSEs of the FFT 
are even better, exhibiting a constant MSE of 0.4 for block 

25 sizes up to 4096 bits. When the demonstrated performance of 
the present invention is viewed in light of the present 
capabilities and increasing need for longer block sizes, the 
advantages of the present invention become clear. 
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5 



N 


IntDCT-IV 


IntlDCT-IV 


IntFFT 


IntlFFT 


8 


0.537 


0.537 


0.456 


0 . 371 


16 


0.546 


0.546 


0.480 


0 . 412 


32 


0.549 


0.548 


0.461 


0 .391 


64 


0.550 


0.550 


0.4 62 


0 . 393 


128 


0.551 


0.551 


0. 461 


0 .391 


256 


0.552 


0.552 


0.461 


0-391 


512 


0.552 


0. 552 


0.461 


0 .391 


1024 


0.552 


0.552 


0. 460 


0 .391 


2048 


0.552 


0.552 


0. 461 


0 .391 


4096 


0.553 


0.552 


0.461 


0-391 



10 



Table I 
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